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Abstract.  We  study  the  scaling  properties  of  heterogeneities  in  nematic  (liq¬ 
uid  crystal)  polymers  that  are  generated  by  pressure-driven,  capillary  Poiseuille 
flow.  These  studies  complement  our  earlier  drag-driven  structure  simulations 
and  analyses.  We  use  the  mesoscopic  Doi-Marrucci-Greco  model,  which  incor¬ 
porates  excluded-volume  interactions  of  the  rod-like  particle  ensemble,  distor- 
tional  elasticity  of  the  dispersion,  and  hydrodynamic  feedback  through  orien¬ 
tation  dependent  viscoelastic  stresses.  The  geometry  likewise  introduces  an¬ 
choring  conditions  on  the  nano-rods  which  touch  the  solid  boundaries.  We 
first  derive  flow-orientation  steady-state  structures  for  three  different  anchor¬ 
ing  conditions,  by  asymptotic  analysis  in  the  limit  of  weak  pressure  gradient. 

These  closed-form  expressions  yield  scaling  laws,  which  predict  how  lengthscales 
of  distortions  in  the  flow  and  orientational  distribution  vary  with  strength  of 
the  excluded  volume  potential,  molecule  geometry,  and  distortional  elasticity 
constants.  Next,  the  asymptotic  structures  are  verified  by  direct  numerical 
simulations,  which  provide  a  high  level  benchmark  on  the  numerical  code  and 
algorithm.  Finally,  we  calculate  the  effective  (thermal  or  electrical)  conductiv¬ 
ity  tensor  of  the  composite  films,  and  determine  scaling  behavior  of  the  effective 
property  enhancements  generated  by  capillary  Poiseuille  flow. 
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1.  Introduction 

The  rheological  behavior  of  Poiseuille  flow  of  liquid  crystal  polymers  is  of  interest 
for  technology  applications  in  processing  of  high  performance  fibers  and  films.  Rey 
and  co-workers  have  studied  capillary  Poiseuille  flows  using  the  Leslie-Ericksen 
theory  for  discotic  nematic  liquid  crystals  ([4,  5,  6]).  Denniston  et  al.  have  used 
lattice  Boltzmann  simulations  of  Landau-deGennes  orientation  tensor  models  to 
explore  the  behavior  of  liquid  crystal  polymers  subject  to  Poiseuille  flow  ([9]).  In 
this  work  we  study  capillary  Poiseuille  flow  using  the  Doi-Marrucci-Greco  model, 
which  also  employs  an  orientation  (second-moment)  tensor  description  of  the  rod 
distribution. 

We  extend  our  early  work  ([3],  [17])  in  the  following  ways:  1)  asymptotic  results 
for  tilted  anchoring  condition  are  derived;  2)  the  rotary  diffusivity  is  explicitly 
coupled  in  the  asymptotic  results;  3)  direct  numerical  simulations  are  carried  out 
to  validate  asymptotic  results;  and  4)  effective  conductivity  properties  of  Poiseuille 
flow-generated  composite  films  are  calculated. 
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The  paper  is  organized  as  follows.  First,  we  describe  the  model  of  liquid  crystal 
polymer  hydrodynamics  proposed  by  Doi,  Marrucci  and  Greco  ([10,  25,  26,  27]). 
Second,  we  present  the  structure  scaling  properties  of  the  orientation  tensor  for 
different  anchoring  conditions;  this  is  achieved  by  asymptotic  analysis  that  yields 
exact  orientation  modes  with  spatial  variations  controlled  by  material  and  boundary 
conditions  in  the  limit  of  weak  pressure  gradient.  We  perform  direct  numerical 
simulations  to  verify  the  asymptotic  results,  and  then  observe  new  properties  as 
the  asymptotic  conditions  are  violated.  Finally,  we  calculate  the  effective  (thermal 
or  electrical)  conductivity  tensor  of  the  composite  films,  and  determine  scaling 
behavior  of  the  effective  property  enhancements  generated  by  capillary  Poiseuille 
flow. 


2.  Model  Formulation 

We  consider  capillary  Poiseuille  flow  of  nematic  liquid  crystal  polymers  (LCPs) . 
Capillary  Poiseuille  flow  is  the  flow  between  two  non-slip  boundaries  ai  y  =  —h 
and  y  =  h  driven  by  a  pressure  gradient  in  this  context.  The  flow  is  described 
by  a  velocity  field  v  =  {vj;{y,t),0,0)  with  the  centerline  of  the  pipe  located  at 
y  =  0.  Figure  1  depicts  the  cross-section  of  the  Poiseuille  flow  on  the  (x,  y)  plane. 
For  simplicity  here  we  suppress  variations  in  the  direction  of  flow  (x)  and  primary 
vorticity  direction  (z),  and  transport  in  the  vertical  (y)  direction. 


Figure  1.  Plane  Poiseuille  flow  geometry.  Non-slip  boundary  con¬ 
ditions  for  the  velocity  and  boundary  anchoring  for  the  orientation 
tensor  given  by  the  stable  nematic  rest  state  are  prescribed,  with 
major  director  angle  i/>o  =  0  shown  here. 


The  dimensionless  governing  equations  consist  of  the  balance  of  linear  momen¬ 
tum,  stress  constitutive  equation,  continuity  equation,  and  the  equation  for  the 
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orientation  tensor  ([17]): 

^  =  V-i-pl  +  r),  p  =  ex, 

=  (ifs +M3)D  +  aa7^(Q) 

+  {AQ  :  Q(Q  +  f)  -  i(AQQ  +  QAQ)  -  lAQ} 

+  3f^  {i(QAQ  -  AQQ)  -  i(VQ  :  VQ  -  VVQ  :  Q)} 


(1) 


+/ii 


(Q+3)D 


r)(Q  +  -) 


M2D  :  Q  (  Q  +  - 


V  •  V  =  0, 

=  llQ-Qll  +  a[DQ  +  QD]  +  fD-2aD:Q(Q+f) 


“i  {^(Q)  + 

where  v  is  the  velocity  field  for  the  flowing  LCP,  p  is  the  pressure,  I  is  the  unit  tensor, 
T  is  the  total  stress,  a  is  a  dimensionless  parameter  depending  on  the  molecular 
aspect  ratio  r  of  spheroidal  molecules 

—  1 


Re  is  the  solvent  Reynolds  number,  Er  is  the  Ericksen  number  which  is  a  reciprocal 
elasticity  constant,  pi  (i  =  1,2,3)  are  three  nematic  Reynolds  numbers,  a  is  a 
normalized  entropic  parameter,  D  is  the  rate-of-strain  tensor  and  H  is  the  vorticity 
tensor, 

(2)  D  =  ^(Vv  +  Vv^),  n=  i(Vv-Vv^). 

Here  the  superscript  T  denotes  the  transpose  of  a  second  order  tensor. 

(3)  E(Q)  =  (1  -  iV/3)Q  -  AQ2  +  AQ  :  Q(Q  +  1/3), 


(4) 


1, 


with  constant  rotary  diffusivity. 


A  = 


(1  —  |Q  :  Q)^,  with  orientation-dependent  rotary  diffusivity. 


The  no  slip  condition  at  the  bounding  surfaces  is  used  for  the  scaled  axial  velocity 


(5) 


Tally— ±1  —  0- 


Homogeneous  mesophase  anchoring  is  assumed  at  the  boundaries,  which  assumes 
that  the  rod  ensemble  is  in  an  equilibrium  bulk  phase.  We  work  at  a  concentration 
where  the  rest  state  is  a  stable  ordered  equilibrium,  called  the  nematic  phase: 

Q|y=±i  =  so(nn-  1), 


(6) 


So  —  4(1  +  aiv  )> 


where  sq  is  the  stable  uniaxial  order  parameter  specified  by  the  nematic  concen¬ 
tration  A  >  |,  and  n  is  the  equilibrium  uniaxial  director.  We  will  fix  A  >  3  so 
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that  the  nematic  phase  is  the  unique  stable  equilibrium,  limiting  the  chances  of 
potential  hysteresis  and  bi-stability  in  flow-induced  structures.  We  assume  n  lies 
in  the  how-flow  gradient  plane  at  some  experimentally  dictated  anchoring  angle  'ipQ 
with  respect  to  the  how  direction, 

(7)  n  =  (cos  ipo,  sin i/'o,  0). 

ipo  =  0,  ipo  =  and  0  <  ipo  <  ^  correspond  to  tangential,  homeotropic  and  tilted 
anchoring,  respectively. 

We  consider  the  in-plane  mesophase  orientation  of  LCPs  with  two  directors  of  Q 
conhned  to  the  plane  (x,  y),  but  still  admitting  biaxiality.  This  restriction  is  not  eas¬ 
ily  lifted  for  the  analysis  to  follow;  numerically,  there  is  not  limitation.  The  point, 
however,  is  to  glean  some  prediction  of  the  types  of  spatial  structures  that  arise, 
out  of  equilibrium  phases,  due  to  the  combined  inhuences  of  a  pressure  gradient 
induced  how,  and  anchoring  of  the  rod  ensemble  at  the  walls.  We  will  hnd,  consis¬ 
tent  with  our  earlier  analyses,  two  types  of  structure  modes:  extended  structures  in 
the  orientational  distribution  which  accompany  the  Poiseuille  how  geometry;  and, 
boundary  layers  which  arise  due  to  wall  conhnement.  Drilling  somewhat  deeper 
in  these  analytical,  though  asymptotic,  insights  into  structures,  a  more  important 
question  is  a  prediction  of  the  scaling  laws  of  each  type  of  structure  mode.  Namely, 
it  is  impossible  to  explore  this  multi-parameter  space  even  with  large-scale  com¬ 
putations;  guidance  as  to  how  the  lengthscales  of  spatial  heterogeneity  vary  with 
parameters  is  of  paramount  value  in  guiding  our  subsequent  numerical  simulations. 
Thus,  we  make  no  apologies  either  for  the  reduction  in  physical  space  dimensions 
from  three  to  one,  or  for  the  reduction  in  orientational  degrees  of  freedom,  namely 
from  a  fnll  probability  distribution  function  with  inhnite  degrees  of  freedom  to  a 
second-moment  tensor  with  generally  hve  degrees  of  freedom,  and  now  to  a  second- 
moment  whose  peak  axis  of  orientation  must  lie  in  the  plane  of  the  primary  flow 
and  flow-gradient  directions.  These  assumptions  shall  lead  to  explicit  structure 
modes  and  scaling  laws,  in  flow  and  particle  orientation,  as  well  as  an  indication  of 
when  these  asymptotic  approximations  will  break  down.  Those  insights  will  then 
indicate  grid-scale  resolution  of  algorithms  across  parameter  space. 

The  in-plane  orientational  constraint  implies  that  two  components  in  the  Carte¬ 
sian  representation  of  Q  must  vanish,  Q^z  =  Qyz  =  0.  Alternatively,  the  orientation 
tensor  can  be  written  in  terms  of  directors  and  order  parameters: 

(8)  Q  =  s{y,  t)(nn  -  ^)  +  /?(j/,  t){n^n-^  -  ^ ), 

with  the  directors  n,  confined  to  the  {x,  y)  plane  and  parametrized  by  the  in¬ 
plane  Leslie  angle  'ip{y,t), 

(9)  n  =  (cos  V',  sin^,  0),  =  (—  sm-i/j,  costp,  0). 

The  third  director  is  rigidly  constrained  along  the  vorticity  axis. 

The  explicit  coordinate  change  between  these  Q  representations  is: 

Qxx  =  s(cos^  /3(sin^  ^  ~  ’ 

Qxy  =  {s  —  (3)  sin  Ip  cos  ip  , 

Qyy  =  s{sin^ip-  ^)  +  (3{cos^  ip  -  . 


(10) 
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With  the  biaxial  representation  (8)  of  Q,  the  tensor  equation  in  (5)  can  be 
written  explicitly  in  terms  of  the  two  order  parameters  (s,  j3)  and  the  Leslie  angle 


ds 


1 


2N 


dip  ^ 


1 

9Er 


3“^  s  3“^  3 

-(l-s)(l-/3  +  2s)^  +  s(l-s  +  2/3)^ 


,  9Voo  ,  , 

+  0  -^9{s,P)  sm2V>, 
3  dy 


dl3 


1 


(11) 


2N 


dip. 


^  =  -X  j  UW  +  —sPil  +  s-P)  +  —iP-s)  gid,  s)  (-) 


9Er' 


1 

9Er 


«1-<3  +  2»)0-(1-4)(1-»  +  2/3)0 


dip 

dt  6(s  —  /3)  \  dy 


1  r  ^1/ 

<  {3f3  —  3s  +  a(2  +  s  +  /?)  cos  2ip) 


(«  dy\ 


where 

(12) 


C/(s)  =  s(l-f(l-s)(2s  +  l)), 
g{s,  /3)  =  1  +  3s/3  —  /3  +  2s  —  3s^  . 


The  momentum  equation  is  reduced  to  a  single  equation  for  the  axial  velocity 
component  Vx, 


dvx  _  dTxy 
dt  dy 


-  e, 


C/(s)  -  C/(/3) 


4iVs/3 


(/3-s) 


sin  2ip 


aa 

ISEr 


(13) 


dip. 


-Ks,d)  TTw  +  h{f^,s)Tr^  +  2(s  -  /3)7(s,/3)(  — ) 


ay- 


ay 


dy' 


sin  2ip 


y  (s  +  /3  +  2)  +  y  (s  -  a)^(l  -  COS#)  +  ^  +  Y 


dvx 

dy 
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where 


(14) 


=  g{s,(3)  +g{f3,s), 
h{s,  P)  =  {1  -  P  +  2s){l  +  (3  —  s). 


We  seek  steady  solutions  of  the  dimensionless  governing  equations  (11-14)  subject 
to  the  boundary  conditions 


(15)  Vx\y=:izl  0?  '^ly— ±1  P\y—±1  O5  'P\y=±l  V^O- 

We  posit  a  formal  asymptotic  expansion  in  the  small  e  limit,  consistent  with  the 
above  boundary  conditions: 


(16) 


s  =  ^sW(^/)e^  P  =  Y,P^>^\y)eK 

k=0  k=l 

Alternatively,  the  Cartesian  representation  of  Q  is  expanded  in  the  form: 

00 

(17)  Q  =  ^Q(")(y)e", 

fc=0 


with  a  quiescent,  homogeneous  equilibrium  given  by  (6).  and 

the  components  of  are  explicitly  related  by  : 

_  2  Q^Jjj  cos  2po  -  (Q^y  -  Q^yJ )  sin  2po 

2^  ’ 

=  ^  (Qix  +Qm)  +  \  (Qxx  -  QyJ )  cos  2iPo  +  g  W  sin  2Po , 

=  ^  {Qxx  +Qm)-\  (Qxx  -  QyJ )  cos  2Po  -  giV  sin  2i/>o  • 

We  remark  that  anchoring  conditions  at  the  two  boundaries  are  assumed  iden¬ 
tical.  Different  conditions  at  each  boundary  require  a  non-homogeneous  steady 
structure  at  leading  order  ([34]),  which  we  will  not  pursue  for  this  study. 


(18) 


3.  Steady  Asymptotic  Structures 

From  either  representation  of  Q,  the  above  expansions  are  inserted  into  the  full 
system  of  flow-nematic  equations,  (l)-(4)  or  (11)-(12),  equations  for  the  0(e)  and 
O(e^)  variables  derived,  and  then  systematically  solved.  We  now  state  the  results, 
following  the  notation  in  expansion  (16). 


3.1.  Velocity  Structure.  The  axial  velocity  at  0(e),  vi^\  reproduces  the  simple 
Poiseuille  flow  structure'. 


(19) 


(1  -  y^), 
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where 


^  +  /ii(2  +  So)  +  3/43  +  |/42Sq  sin^(2'i/;o)  +  (Mi M2  +  M^)a 

3A 

Ml  =  — - (1  -  Al  cos2V’o), 

2  +  So 

M2  =  3so  (1  -  Ai  cos27/>o), 

M3  =  sin^(2i/'o)  A  (2  +  so  -  3so), 

a(2  +  So) 


3.2.  Director  Structure.  The  first  consequence  of  these  exact  mesoscopic  so¬ 
lutions  is  quite  striking:  the  director  is  seen  to  only  participate  in  the  long-range 
extended  structure  and  to  not  vary  in  the  boundary  layers.  From  equation  (18),  the 
director  angle  ^|)  is  explicitly  given  by 

(21)  Ip  =  tpo eCpp  {1  —  y'^) Cp  =  CyMiEr. 

By  inference,  the  confinement-induced  boundary  layer  structures  are  therefore  re¬ 
sponsible  only  for  focusing  of  the  orientational  distribution. 

It  is  easy  to  show  that  Mi  M2  -I-  M3  >  0.  Then,  in  the  denominator  of  C^,  all 
terms  are  positive  definite,  so  that  is  positive  definite.  Mi,  however,  is  not: 
it  depends  on  whether  the  Leslie  parameter  Al  is  greater  or  less  than  one.  As  in 
Couette  (plate-driven)  flow  ([17]),  this  controls  the  direction  of  twist  of  the  major 
director  from  plate  to  plate,  from  formula  (21). 

3.3.  Order  Parameter  Structure.  We  now  extract  degree  and  scales  of  molec¬ 
ular  elasticity  induced  by  slow  Poiseuille  flow.  The  striking  feature  is  extreme 
sensitivity  to  anchoring  angle.  From  the  s,  (3  formulas,  (18),  it  follows  that  and 

vanish  if  and  only  if  -i/jq  =  0,  |  (mod  tt).  These  were  the  sole  boundary  condi¬ 
tions  considered  in  ([3]),  so  the  tilted  structures  presented  here  will  have  different 
features,  again  matching  the  similar  phenomenon  discovered  in  Couette  flow  ([17]). 

•  Tangential  anchoring  (gpQ  =  0) 

Order  parameter  distortions  enter  at  O(e^)  through  two  boundary  layers  and  the 
long-range  permeation  structure: 

s  =  so-ke^s^^^-kO(e^), 

/3  =  e2/3(2) 


/3(2)  =  c[^)  cosh{Er^/‘^Gy)  +  -k  cf  ^(2)^ 

sO)  =  2Cf'^  cosh{Er^/'^Hy)  +  2Cf'^  cos\v{Er^/‘^G  y) 


(22) 
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where 


4(2)  2C'p(l  —  so)(3soC'p  —  2aCvErAcos2'tl;o) 

2  ^  Er{3  -N  +  2Nsl  +  2Nso)  ’ 

^(2)  12(^2^^  —  CpSQ  +  aCyCpEr  Acos2'ipo) 

^2  =  -^2  - 

(2)  ^  2(cf )  +  C^pSp)  (2)  ^  )  +  e!^^ 

2  ErG^  ’  1  cosh(£;ri/2G)  ’ 

(2)  ^  soiG^+6N)c[^'> 

^  (l+2so)(G2-i?2)’ 

(2)  _  6[ErNsoA^^'>  +  Cp  (3so  +  1)  (3soGp  -  2aG^ErA  cos  2i/>o)] 

^  “  ErH^{l  +  2so) 

+  +  EtNG^^^) 

+2Gp  (1  +  3so)  (GpSo  —  aGt,i?rAcos2i/)o)] 

(2)  2(1  +  2so)cf  ^  -  2so(2Gf  ^  +  SErNE^^^)  -  2G^soil  +  3so) 

“  ErH‘^{l  +  2so) 

(2)  _  2gP  cosh{ErE^G)  +  ^  +  cjf 

^  2  cosh(i?ri/2/f) 

Combined  with  the  director  results,  (21),  we  conclude  that  tangential  anchoring 
leads  to  a  director-dominated,  non-uniform,  long-range  structure  with  Er~^  mean 
scaling,  the  same  behavior  as  in  Couette  flow.  This  implies  that  pressure-driven 
and  drag-driven  flows,  which  are  respectively  parabolic  and  linear  shear  layers,  do 
not  alter  the  molecular  elasticity-dominated  orientational  gradients. 

•  Normal  (homeotropic)  anchoring  (ipo  =  7r/2) 


The  structure  model  (11)-(15)  for  in-plane  Q  tensors  admits  a  symmetry:  (a,  ip) 
(—a,  I  Eip).  This  property  implies  the  asymptotic  solution  with  ipo  =  77/2  can  be 
obtained  directly  from  (22).  Thus  normal  and  tangential  anchoring  structures  are 
guantitatively  different,  with  changes  in  numerical  pre-factors,  yet  have  gualitatively 
similar  scaling  laws.  The  figures  will  illustrate  these  features. 

•  Tilted  anchoring  if)  <  ipo  <  it  12  or  7r/2  <  <  tt) 

Any  pre-tilt  induces  an  order  parameter  response  of  the  same  amplitude,  0(e), 
as  the  director  distortion,  yet  localized  in  the  two  boundary  layers: 

s  =  So  +  +  o(e^). 


(23) 


/3  =  £/?(') +0(6^), 
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where 

=  2Ai  sinh(Er^'^^Dy)  +  A2  y, 

=2  Cl  si-nh{Er^/^Fy)  -  Ca  y  +  2  Cg  si-nh{Er^/^ Dy), 


B  =  9  +  6Nso-3N  +  6Nsl,  D  = 
6aCvAsin{2ipo)  (1  —  sq) 


B 


1-so 


Ao,  = 


B 


Ai  =  - 


Ao 


2  smh(Er^CE>)  ’ 


Bi  =  ISNsl  —  3N  —  6Nso  +9,  -82  =  6oC^Asin(2i/)o)  (1  +  2so  —  3sg), 


B.  = 


E  = 


QNso{l  —  80)^2  +  B2 
'Bi  ^ 


Bi 


1  +  So  ~  2s, 


2  ’ 
0 


C2  = 


6-/Vso(l  —  so)A2  +  B2 
Hi  ^ 


_  (so  -  Sq)  Ai  {Er  D'^  +  6Er  N) 

^  ~  (l  +  so-2s§-Si  Er)ErD‘^  ’ 

_  ^3  -  2C3  siiAi(Er^/‘^D) 

^  2s\h)i(ECCe) 

This  result  for  Poiseuille  flow  is  qualitatively  similar  to  Couette  flow  in  the  slow 
flow  limit  ([3],  [17]).  It  suggests  a  flow- nematic  structure  mechanism:  molecular 
elasticity  is  amplihed  in  local  boundary  layers  when  the  major  director  is  pinned 
at  an  angle  tilted  with  respect  to  the  flow-flow  gradient  axes.  This  result  forebodes 
a  phenomenon  that  only  occurs  at  much  higher  driving  conditions  or  with  much 
’’softer”  elasticity  constants,  namely  the  appearance  of  defects  in  the  orientational 
distribution.  These  structures  correspond  to  local  regions  of  isotropy  (randomness) 
of  the  orientational  distribution.  In  such  regions,  the  local  principal  axes  of  orienta¬ 
tion  are  not  able  to  align  with  the  primary  flow  axes,  which  causes  sharp  boundary 
layers  with  even  smaller  lengthscales,  on  the  order  of  Er~^  according  to  the  predic¬ 
tions  above.  Since  Er  for  nematic  polymers  are  on  the  order  of  10®,  this  suggests  a 
very  stringent  grid-scale  will  be  required  to  resolve  these  structures,  and,  an  adap¬ 
tive  mesh  refinement  technique  will  be  required  since  there  is  no  a  priori  knowledge 
of  where  these  defects  will  arise.  Finally,  in  nematic  polymers  as  opposed  to  liquid 
crystals,  defects  are  transient  as  opposed  to  preserved,  and  attractors  are  dynamic 
rather  than  predominantly  steady;  this  means  under-resolution  in  time  integrators 
or  simply  a  focus  on  steady  structures  are  essentially  useless. 

4.  Effective  Conductivity  Tensors 

We  now  take  the  results  from  above  on  pressure-driven  processing  flows,  and 
predict  what  kinds  of  conductivity  features  are  associated  with  these  orientational 
distributions.  This  is  part  of  an  effort  within  our  larger  research  group  to  link 
processing  of  nematic  polymer  nano-composites  (NPNCs)  to  effective  properties. 
In  ([32]),  we  developed  a  connection  between  the  probability  distribution  function 
of  nematic  polymers  and  the  effective  conductivity  tensor  of  low-volume-fraction 
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nanocomposite  materials  based  on  homogenization  theory.  One  simply  needs  to 
posit  the  properties  of  the  nano-rods  and  of  the  matrix,  where  each  phase  is  assumed 
to  be  an  isotropic  material.  Extensions  to  anisotropic  nano  and  matrix  phases  are 
straightforward,  but  a  detail  that  only  clouds  the  essential  predictions  we  wish  to 
convey  here. 

This  approach  considers  spheroidal  molecules  (rods  in  this  instance)  with  molec¬ 
ular  orientation  m,  aspect  ratio  r,  isotropic  conductivity  a2  and  volume  fraction 
02  (assumed  quite  small,  indeed  in  practical  nano-composites  on  the  order  of  1%) 
populated  in  a  matrix  of  conductivity  cti  according  to  the  orientation  distribution 
function  /(m)  of  the  ensemble.  Conductivity  obeys  elliptic  equations  which  are  vir¬ 
tually  identical  for  thermal,  electric,  and  dielectric  properties.  The  volume  fraction, 
02,  of  the  nematic  polymer  inclusions  is  proportional  to  the  dimensionless  strength, 
N,  of  the  Maier-Saupe  excluded  volume  potential: 


(24) 


02  = 


Ntt 

”87’ 


where  r  =  l/d  is  the  aspect  ratio  of  the  molecular  spheroids  with  length  I  and 
diameter  d.  Then,  the  effective  conductivity  tensor  in  closed  form  is  derived  in  the 
low- volume- fraction  limit  {02  <<  1)  as: 


Cr®  =  CTiI  -I-  ai02{(72  —  (Ji) 


(25) 


-f 


fL2  +  CTi  -  ((72  -  0-i)La 
((72  -  0-i)(l  -  3La) 


[((72  +  CTi)  -  ((72  -  ai)La][ai  +  (((72  -  (7i)La] 


M(/)  +  Oi0i), 


where  I  is  the  3  by  3  identity  matrix,  La  is  the  spheroidal  depolarization  factor 
depending  on  the  nematogen  aspect  ratio  r  through  the  relation 


(26) 


La  = 


1-e" 


—  In 
2e 


1 


1  -  e 


-  1 


=  a/I  -  r-2. 


Here  the  orientation  probability  distribution  function  of  the  LCP  inclusions,  /,  is 
governed  by  the  Smoluchowski  equation  of  the  Doi-Hess  kinetic  theory  for  quiescent 
or  flowing  nematic  polymers  ([18]);  M(/)  is  the  second-moment  of  the  orientation 
probability  density  /.  Recall  that  the  orientation  tensor  Q  of  the  previous  devel¬ 
opment  is  simply  Q  =  M  —  ll. 

From  equation  (25)  it  is  obvious  that  the  effective  conductivity  shares  the  same 
three  principal  axes  with  the  second  moment  tensor  M(/).  We  denote  the  three 
principle  values  of  conductivity  (i.e.  eigenvalues  of  ct®)  as  cr®  >  (7|  >  (7|,  where 
the  largest  principal  value,  cr®,  is  the  maximum  effective  conductivity.  The  relative 
principal  value  enhancements  are  given  by 

e  _  _ 

(27)  ^ i  =  1,2,3, 

where  £1  >  £2  >  £3-  The  value  £1  —  £2  depicts  property  anisotropy  with  respect  to 
the  two  principal  conductivity  directions,  whereas  £1  captures  the  maximum  scalar 
enhancement  due  to  the  nematic  polymer  inclusions.  For  example,  £1  =  1  implies 
cr®  =  2cti  or  100%  gain  of  the  bulk  composite  conductivity. 

Heterogeneous  effective  properties  of  sheared  and  extension-enhanced  nematic 
polymer  nanocomposites  have  been  studied  in  ([18,  33]).  In  the  next  section  we 
will  study  the  anisotropy  and  effective  properties  enhanced  by  Poiseuille  flows. 
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5.  Numerical  Results 

5.1.  Numerical  Methods.  In  order  to  numerically  find  the  steady  state  solutions 
of  the  boundary  value  problem  (11-15),  we  use  Newton’s  method  to  iteratively  ob¬ 
tain  the  steady  state  solutions.  Our  codes  cannot  resolve  parameter  regimes  where 
Er  is  above  a  critical  threshold  value;  this  may  be  either  due  to  singular  structures 
or  due  to  structure  transitions  where  steady  structures  give  rise  to  transient  at¬ 
tractors.  In  the  high  Ericksen  number  limit,  where  the  lengthscales  of  distortional 
elasticity  are  predicted  to  scale  like  Er~^  and  Er~^l’^ ^  adaptive  gridding  will  be 
required  since  one  does  not  know  a  priori  where  these  internal  boundary  layers  will 
form.  Some  of  these  numerical  tools  have  been  developed  in  our  group,  but  none 
have  yet  been  applied  to  pressure-driven  flows. 

5.2.  Numerical  validation  of  the  asymptotic  results  and  vice  versa.  Fig¬ 
ures  2-4  compare  the  exact  asymptotic  formulas  for  the  director  angle  order 
parameters  s,  /?,  and  flow  velocity  Vx,  with  direct  numerical  solutions  of  the  steady, 
boundary  value  problem.  Figure  4  also  gives  shear  stress  and  light  scattering  in¬ 
tensity  predictions,  which  provide  a  basis  for  experimental  validation.  For  the 
computations  shown  here,  we  fix  a  concentration  in  the  nematic  regime  N  =  6  and 
the  aspect  ratio  r  =  3  (i.e.,  a  =  0.8),  consistent  with  our  earlier  studies  ([14,  17]). 

We  begin  with  confirmation  of  the  formulas  for  0  <  e  ^  1  and  0  <  e  •  Er  <C 
1,  Figure  2,  where  the  asymptotic  formulas  are  nearly  identical  to  the  numerical 
solutions.  For  tangential  and  normal  anchoring,  we  set  e  =  0.1,  Er  =  1;  for  tilted 
anchoring,  '00  =  f  j  we  have  to  lower  e  =  0.05  with  the  same  Er  =  1  to  maintain 
agreement  between  the  asymptotic  formulas  and  the  numerical  solutions. 

Next  we  illustrate  breakdown  of  the  asymptotic  formulas  with  direct  numerical 
solutions  when  either  of  the  two  asymptotic  conditions,  e  1  and  e  ■  Er  1,  are 
gradually  pushed  out  of  the  asymptotic  regime.  Our  purpose  here  is  to  monitor 
emergence  of  new,  non-asymptotic  structures,  and  to  identify  sources  of  additional 
scaling  behavior. 

Figure  3  imposes  tangential  (Figure  3c  and  d)  and  normal  anchoring  (Figure  3e 
and  f)  with  e  =  0.3,  Er  =  10.  For  tilted  anchoring.  Figures  3a,  3b,  a  more  stringent 
e  •  Er  condition  is  needed  for  asymptotic  agreement,  with  e  =  0.5,  Er  =  1.  One 
finds  quantitative  asymptotic  accuracy. 

In  Figure  4a,b,  we  raise  the  pressure  gradient  to  (e  =  0.5),  with  Er  =  20  (for 
which  e  •  Er  =  10),  and  present  results  only  with  normal  anchoring.  Note  that  the 
velocity  profile  departs  from  being  parabolic  near  the  boundaries.  The  concavity 
of  the  velocity  profile  changes  from  concave-up  to  concave-down  and  then  concave- 
up  again.  We  further  give  predictions  of  shear-gap  stresses  (Figure  4c)  and  light 
intensity  patterns  (Figure  4d)  associated  with  these  steady  structures.  Neither  Ni 
nor  N2  changes  sign  across  the  rectangular  pipe.  The  light  scattering  intensity 
across  the  pipe  associated  with  the  numerical  solutions  is  given  in  Figure  4d.  This 
intensity  pattern  can  be  viewed  as  a  combination  of  pure  nematic  patterns  of  the 
Doi-Marrucci-Greco  theory,  ([15,  16]),  where  the  flow  coupling  has  led  to  selection 
of  order-parameter  structures  near  the  plates  and  director  structures  in  the  interior. 

5.3.  Effective  conductivity  properties.  Our  main  focus  in  this  subsection  is 
the  determination  of  effective  properties  associated  with  capillary  Poiseuille  flow. 
We  consider  typical  model  NPNCs  with  the  nanorod  aspect  ratio  r  =  100  (and 
one  with  r  =  3  for  comparison)  and  the  conductivity  contrast  between  the  matrix 
solvent  and  nanophase  as  cti/ct2  =  10“^.  This  ratio  is,  in  fact,  quite  modest; 
some  composites  have  contrasts  of  ten  or  greater  orders  of  magnitude.  Figure  5 
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Figure  2.  Comparison  of  the  numerical  solutions  (circles  or  trian¬ 
gles)  of  the  steady  boundary  value  problem  and  the  exact  asymp¬ 
totic  formulas  (solid  or  dashed  lines  ).  Top  row:  tilted  anchoring 
{ipo  =  f )  with  e  =  0.05  and  Er  =  1.  Middle  row:  parallel  an¬ 
choring  (ipQ  =  0)  with  e  =  0.1  and  Er  =  1.  Bottom  row:  normal 
anchoring  {t/jo  =  f )  with  e  =  0.1  and  Er  =  1.  In  all  figures,  circles 
and  solid  lines  use  the  left  (vertcal)  axis;  diamonds  and  dashed 
lines  use  the  right  axis. 


shows  the  maximum  scalar  conductivity  enhancement  (cmax)  across  the  gap  for 
three  different  anchoring  conditions.  The  top  row  corresponds  to  the  aspect  ratio 
r  =  100  whereas  for  the  bottom  row  r  =  3.  It  is  apparent  that  one  can  achieve 
different  effective  conductivity  enhancements  across  the  gap  by  tuning  anchoring 
conditions.  The  nematogen  aspect  ratio  effect  is  also  significant  from  this  figure. 
The  larger  the  aspect  ratio,  the  more  enhancement  in  effective  conductivity. 

It  is  also  interesting  to  note  that  the  shape  of  the  solution  for  tangential  anchoring 
is  sensitive  to  the  value  of  aspect  ratio.  In  Figure  6  we  plot  the  maximum  effective 
conductivity  enhancement  and  anisotropy  versus  Poiseuille  flow  rate  (Figure  6a  and 
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X  10"' 


(e) 


(f) 


Figure  3.  Comparison  of  the  numerical  solutions  (circles  or  trian¬ 
gles)  of  the  steady  boundary  value  problem  and  the  exact  asymp¬ 
totic  formulas  (solid  or  dashed  lines).  Top  row:  tilted  anchoring 
{ipo  =  f )  with  e  =  0.5  and  Er  =  1.  Middle  row:  parallel  anchor¬ 
ing  ('00  =  0)  with  e  =  0.3  and  Er  =  10.  Bottom  row:  normal 
anchoring  (0o  =  f )  with  e  =  0.3  and  Er  =  10. 

c)  and  volume  fraction  (Figure  6b  and  d)  for  three  different  anchoring  conditions. 
The  parameters  used  here  are  r  =  100  and  Er  =  1.  In  Figure  6a  and  c,  the  volume 
fraction  is  fixed  to  be  02  =  0.0236  and  we  vary  the  Poiseuille  flow  rate.  In  this 
case,  tilted  anchoring  gives  the  largest  conductivity  enhancement  and  anisotropy. 
In  Figure  6b  and  d,  the  flow  rate  e  is  fixed  to  be  0.1  and  we  vary  the  volume 
fraction  62-  It  turns  out  that  the  conductivity  enhancement  and  anisotropy  are 
both  increasing  functions  of  62  and  insensitive  to  anchoring  conditions. 

6.  Conclusions 

We  have  studied  the  structure  scaling  properties  of  nematic  liquid  crystal  poly¬ 
mers  under  pressure-driven  Poiseuille  flow  both  analytically  and  numerically.  A 
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(a) 


(b) 


(c) 


Figure  4.  Top  row:  Comparison  of  the  numerical  solutions  (cir¬ 
cles  or  diamonds)  of  the  steady  boundary  value  problem  and  the 
exact  asymptotic  formulas  (solid  or  dashed  lines)  for  normal  an¬ 
choring  (i/>o  =  7r/2)  with  e  =  0.5  and  Er  =  20.  Bottom  row: 
Rheological  structure  profiles  associated  with  the  numerical  re¬ 
sults  of  the  top  row.  Here  we  provide  the  first  (A^i  =  th  —  T22) 
and  second  {N2  =  T22  —  T33)  normal  stress  differences,  the  shear 
stress  (ri2),  and  the  light  scattering  intensity  function  ([28]), 
I  =  Jo  sin^  (27101  s),  across  the  plate  gap,  where  Oi  is  the  ratio 
of  the  sample  thickness  to  the  wavelength  of  incident  light.  The 
normalized  birefringence  s  measures  the  difference  between  orien¬ 
tation  along  the  flow  direction  and  along  the  vorticity  axis. 


rich  phenomenology  is  apparent  because  of  the  coupling  between  the  rod  ensemble 
(through  excluded  volume  and  distortional  elasticity)  and  the  flow  field.  When  the 
normalized  pressure  gradient  e  and  the  scaled  Ericksen  number  e-Er  are  both  small, 
the  mesoscopic  Doi-Marrucci-Greco  model  predicts:  an  Er~^/^  scaling  associated 
with  boundary  layer  modes  residing  only  in  the  order  parameters;  and,  an  Er~^ 
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Figure  5.  Effective  conductivity  enhancements  for  3  different  an¬ 
choring  conditions.  Top  row:  r  =  100;  Bottom  row:  r  =  3.  First 
column:  parallel  anchoring;  Second  columm:  normal  anchoring; 
Third  column:  titled  anchoring  (i/>o  =  f  )•  Here  e  =  0.5,  Er  =  1. 


mean  scaling  law  across  the  rectangular  pipe  residing  in  the  principal  axes  of  ori¬ 
entation.  The  self-consistent,  leading  order,  flow  velocity  is  a  parabolic  Poiseuille 
flow,  attenuated  by  an  explicit  prefactor  which  accounts  for  the  stored  energy  in 
the  rod  dispersion.  Anchoring  conditions  play  an  important  role  in  the  amplitude 
response  of  order  parameters,  which  describe  the  relative  focusing  or  defocusing  of 
the  rod  ensemble  distribution  function. 

This  paper  presents  a  combined  approach  to  understanding  flow-induced  behav¬ 
ior  of  rod-like  macromolecular  dispersions,  namely  blending  modeling,  asymptotic 
analysis,  homogenization  theory,  numerical  analysis,  and  scientific  computation. 
This  is,  however,  just  the  beginning  in  an  overall  strategy  to  design  nano-composite 
materials  from  the  bottom  up.  The  ideal  strategy  consists  of  a  control  wrapper  code, 
where  one  is  specifying  performance  features  such  as  current  distributions  across 
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Figure  6.  Effective  conductivity  enhancement  and  anisotropy 
versus  Poiseuille  flow  rate  e  (a  and  c)  or  volume  fraction  O2  (b 
and  d)  for  3  different  anchoring  conditions. 


a  material  during  operating  conditions.  Using  composition  information,  process¬ 
ing  controls,  and  effective  property  characterizations,  the  challenge  is  to  be  able 
to  modify  composition  and  processing  controls  to  achieve  a  specified  performance 
target.  It  is  this  challenge  that  we  have  only  begun  to  put  into  a  framework,  and 
where  the  guidance  of  modern  applied  mathematics  can  play  a  critical  role.  This 
paper  is  dedicated  to  Max  Gunzburger,  whose  work  and  consultation  over  the  years 
has  inspired  our  group  to  envision  such  a  control  strategy  in  the  first  place,  and 
whose  contributions  will  be  a  key  to  our  anticipated  success. 
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